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SKYWAVE  RADAR  DETECTABILITY  OF  VOLCANIC  AEROSOLS 


S.J.  Anderson 


SUMMARY 

During  1982  there  occurred  several  incidents  in  which  com¬ 
mercial  jet  aircraft  suffered  engine  failure  over  Indonesia 
as  they  passed  through  clouds  of  dust  injected  into  the 
stratosphere  by  an  erupting  volcano  on  Java.  This  paper 
presents  some  theoretical  estimates  of  the  detectability  of 
such  clouds  using ,skyvave  radar.  The  results  imply  that 
detection  would  not  be  possible  with  the  JINDALEE  Stage  B 
radar. 
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1.  INTRODUCTION 

On  24  June  1982  a  Boeing  747  operated  by  British  Airways  suffered  engine  failure 
at  35  000  ft  over  Indonesia  after  flying  through  volcanic  dust  injected  into  the 
atmosphere  by  Mt  Galunggung,  a  7500  ft  volcano  located  in  SW  Java,  which  had  been 
erupting  intermittently  since  5  April  1982.  Subsequently,  on  14 'July  1982,  a 
Singapore  Airlines  Boeing  747  experienced  failure  of  two  engines  and  overheating 
of  a  third  when  it  encountered  dust  at  39  000  ft.  The  air  route  concerned  lies 
partly  within  the  coverage  of  the  JINDALEE  Stage  B  radar  and  the  question  arises: 
could  the  JINDALEE  radar  detect  or  infer  the  presence  of  volcanic  dust  clouds  so 
as  to  alert  air  traffic  to  possible  hazards? 

There  are  several  ways  in  which  volcanic  eruptions  could  conceivably  modify  the 
physical  environment  as  observed  by  a  skywave  radar.  The  most  apparent  of  these 
are: 

(i)  injection  of  volcanic  material  into  the  upper  atmosphere  where  it 
could  reflect  skywave  radar  signals  via  the  Rayleigh  scattering  mechanism; 

(ii)  physical  and  chemical  processes  involving  volcanic  dust  whereby  the 
normal  ionic  equilibria  in  the  upper  atmosphere  are  upset  to  the  point  where 
the  refractive  index  is  appreciably  altered; 

(iii)  generation  of  acoustic-gravity  waves  which  could  propagate  to  iono¬ 
spheric  heights  where  they  would  manifest  themselves  as  travelling  ionospheric 
disturbances. 

Of  these  only  (i)  and  (ii)  are  relevant  to  the  detection  of  dust  clouds. 

This  note  presents  some  rough  estimates  of  the  magnitudes  and  hence  the  signific¬ 
ance  of  these  mechanisms. 


2.  CHARACTERISTICS  OF  VOLCANIC  AEROSOLS 


2.1  Composition 

Individual  volcanoes  vary  significantly  in  the  composition  of  their  ejecta  but 
the  fine,  particulate  tephra  fall  predominantly  into  four  classes(ref . 1) : 

(a)  lithic  fragments  -  particles  of  the  rock  which  made  up  the  summit  of 
the  volcano  prior  to  eruption; 

(b)  pumice  -  a  rock  froth  formed  by  the  rapid  quenching  of  magma,  composed 
of  volcanic  glass,  crystals  of  several  different  minerals,  and  gas  bubble 
voids; 

(c)  crystals  and  crystal  fragments  of  various  minerals  which  were  derived 
from  lithic  fragments  as  well  as  from  the  crystallising  magma;  and 

(d)  shards  of  volcanic  glass  broken  from  the  frothy  magma. 

The  most  important  volcanic  sources  of  high  altitude  contamination  are  the 
so-called  "Plinian"  eruptions,  that  is,  gas-rich,  sustained  explosive  erup¬ 
tions  with  exceptionally  high  (>  10  km)  eruption  columns  (Table  1).  For 
example,  the  plinian  eruption  column  was  the  principal  means  of  transporting 
particles  and  gases  into  the  stratosphere  during  the  18  May  1980  eruption  of 
Mount  St  Helens  (Tables  2  to  4) . 


The  gases  and  nonsilicate  aerosol  particles  contained  in  the  plinian  column 
have  six  possible  sources: 

(i)  tropospheric  air  entrained  by  the  ascending  column; 

(ii)  geothermal  fluids  near  the  magma  body; 

(iii)  vaporisation  of  ground  water  from  melted  snow,  ice  and  rain; 

(iv)  reaction  of  vegetation  with  a  cloud  of  hot  ash  and  gas; 

(v)  gases  released  explosively  from  the  silicate  liquid  in  the  erupted 
magma;  and 

(vi)  gases  released  explosively  from  the  silicate  liquid  in  the  magma 
that  was  not  erupted. 

2.2  Chemistry 

The  effects  of  the  chemistry  of  volcanic  emissions  on  the  stratosphere  and 
lower  ionosphere  are  strongly  dependent  on  the  emission  fluxes  and  the  result¬ 
ing  transport  processes.  It  appears  that  the  most  important  chemical 
processes  are: 

(i)  the  chemical  production  of  aerosols  which  affect  the  Earth's  radia¬ 
tion  budget;  in  particular  the  sulphur-containing  gases  and  water  vapour 
are  important; 

(ii)  the  chemical  interaction  of  volcanic  gases  (and  perhaps  volcanic 
aerosol  surfaces)  with  the  free-radical  ozone  chemistry  of  the  stratosphere. 
The  important  species  are  H2O  and  SOj  and  perhaps  reactive  chlorine  com¬ 
pounds,  CO,  and  oxides  of  nitrogen. 

Most  of  the  emissions  of  sulphur  compounds  eventually  form  sulphate  aerosols. 
The  increasing  dominance  of  sulphate  particles  over  ash  particles  with  time 
in  the  stratosphere  can  be  explained  by:  (a)  early  scavenging  of  ash  by  the 
condensation  of  emitted  water  vapour  to  form  rapidly  settling  ice  crystals; 

(b)  rapid  settling  of  large  ash-containing  particles;  and  (c)  continuing 
productions  of  sulphate  by  oxidation  of  sulphur  gases. 

2.3  Transport  and  dispersion 

In  the  absence  of  strong,  stable  layers  in  the  troposphere  the  eruption  plume 
encounters  no  significant  buoyant  resistance  to  its  upward  progress  until  it 
passes  about  12  km  in  altitude.  Further  penetration  to  heights  up  to  45  km 
has  been  noted  but  25  km  is  more  typical  (Tables  3  and  4;  figure  1). 

The  subsequent  evolution  of  the  plume  is  obviously  critically  dependent  on 
the  stratospheric  winds. 

Inspection  of  SAGE  satellite,  lidar  and  jet  impaction  data(ref.2  to  4) 
suggests  the  preferential  formation  of  thin  laminar  clouds  due  to  wind  shears 
(figures  2  to  4) . 

The  trajectories  of  the  clouds  at  different ( levels  can  meander,  loop  and 
diverge  in  a  very  complex  fashion  making  prediction  extremely  difficult.  This 
is  evident  in  figures  5  and  6(ref.l).  Attempts  to  model  cloud  evolution  have 
not  met  with  great  success.  For  example,  in  the  case  of  the  Caribbean 
Soufriere  eruptions  in  1979,  trajectories  were  calculated  using  a  geostrophic 


ERL-0301-TR 


wind  model  based  on  NOAA  geopotential  maps,  assuming  that  the  injected 
material  resided  on  the  70  mb  surface(ref .4) .  While  a  qualitative  agreement 
was  found  with  subsequent  plume  positions  over  West  Africa  as  observed  by  the 
SAGE  satellite,  the  arrival  times  did  not  match.  Further,  airborne  lidar 
measurements  of  optical  extinction  disagreed  with  the  computed  trajectory 
(figure  7) .  This  discrepancy  is  variously  attributed  to  the  tendency  of 
global  maps  to  emphasise  large-scale  features  of  the  flow  fields  and  neglect 
local  fluctuations,  and  to  the  fact  that  the  geostrophic  approximation  appears 
to  overestimate  wind  speeds  when  used  closer  to  the  equator  than  about  15  of 
latitude. 


3.  ESTIMATES  OF  AEROSOL  EFFECTS 
3.1  HF  radar  cross  section 

The  scattering  of  electromagnetic  waves  from  a  tenuous  distribution  of  small 
particles  is  an  important  problem  which  has  attracted  interest  since 
Lord  Rayleigh's  original  work  in  1899. 

Two  distinct  frameworks  exist  for  dealing  with  "multiple  scattering"  phenomena: 
analytical  theory  and  transport  theory.  Analytical  theory  begins  with  the 
underlying  differential  equations  describing  the  field  properties,  introduces 
the  scattering  and  absorption  properties  of  the  individual  scatterers  and  then 
expresses  the  quantities  of  physical  interest  in  terms  of  appropriate  integral 
equations.  These  can  then  be  solved  by  iterative  or  diagrammatic  techniques. 
Transport  theory  deals  directly  with  the  flow  of  energy  using  the  radiative 
transfer  equation.  A  respectable  attempt  to  calculate  the  transmission  prop¬ 
erties  of  a  volcanic  aerosol  would  require  one  of  these  approaches  but  time 
does  not  permit  such  an  endeavour  on  this  occasion.  Instead  the  simpler 
problem  of  estimating  the  backscattering  radar  cross  section  is  addressed. 

Owing  to  the  extremely  small  dimensions  of  the  aerosol  particles  it  turns  out 
that  the  result  obtained  here  is  equivalent  to  that  which  would  result  from 
the  cumulative  forward-scatter,  single  backscatter  approximation  of 
De  Wolf (ref. 5) . 

Consider  the  contribution  to  the  backscattered  field  from  one  aerosol  particle, 
assumed  spherical.  It  can  be  shown  (reference  6,  page  510)  that 


n2  -  n  k2 a3  -ikR 


,  r  n‘  -  1  -|  k_j 

o  2  ->  R 

_  n  +  2  _ 


with  n  the  refractive  index,  a  the  radius  and  R  the  distance  to  the  field 
point. 

Defining  the  backscatter  cross  section  o^  by 


4*  R2  I  E  J  2 
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The  only  information  available  about  size  distribution  for  this  eruption  is 
that  the  great  majority  of  aerosol  particles  had  a  <  100  Mm  while  about  1% 
of  the  mass  was  present  in  particles  with  a  <  1  Mm. 

From  figure  9  a  reasonable  lower  limit  for  a  theoretical  model  of  N(a)  is 
0.1  Mm.  Any  errors  in  introduced  by  ignoring  particles  smaller  than  this 

are  likely  to  be  negligible  given  the  a6  weighting  in  equation  (4).  Hence  the 
following  constraints  apply: 


1  Mm  . 100  Mm 

a3N(a)  da  =  10-2  /  a3N(a)  da 


(6) 


0. 1  Mm 


0.1  Mm 


4 irp 
3 


100  Mm 


a3  N (a)  da 


5  x  10  6  g  cm  J  (upper  bound) 


-3 


(7) 


0. 1  Mm 


These  two  equations  suffice  to  specify  a  two  parameter  model  of  N(a). 
Physical  considerations  suggest  a  model  of  the  form 


log  N(a)  -  a  log (a)  =  0 


(8) 


with  a  and  0  to  be  determined  from  equations  (6)  and  (7),  and  this  form  is 
supported  by  the  asymptotic  distribution.  Rewriting  equation  (8)  in  exponen¬ 
tial  form,  substituting  in  equation  (6)  and  integrating  yields 
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(9) 


provided  <t  ¥=  -  4. 

Analytic  solution  of  this  equation  is  nontrivial.  Reasoning  that  the  ratio  of 

the  ranges  of  integration  of  the  two  sides  of  equation  (9)  is  about  10~2 ,  so 
the  mass  density  should  therefore  be  only  a  weak  function  of  a  implies  that  u. 
~  -  3.  Using  perturbation  techniques  leads  to 


a  =  -3.024 


(10) 
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Substituting  in  equation  (7) , 


Inserting  equation  (11)  into  equation  (5)  with  D  =  2  km  and  'J',  the  angle  of 
elevation  set  to  6°  yields 


oT0T  =  2.3  x  10-IS  m2  (12) 

This  is  the  radar  cross  section  for  one  square  metre  of  aerosol  measured 
normal  to  the  propagation  vector.  Converting  to  ground  coordinates  and  taking 
radar  parameters  appropriate  for  Mt  Galunggung  observed  with  the  JINDALEE 
radar,  the  effective  radar  cross  section  per  radar  resolution  cell  becomes 

i 

agf£  =  2.4  x  10"8  m2  (13) 

This  value  is  many  orders  of  magnitude  below  detectability  with  a  skywave 
radar. 

It  is  interesting  to  note  that  with  a  microwave  radar  operating  at  10  GHz, 
the  radar  cross  section  per  unit  volume  of  aerosol  would  be  greater  by  a 
large  factor: 


< (microwave)  ~  6.25  x  10* 0  o  (HF) 


(14) 


If  the  beamwidth  of  the  microwave  radar  were  taken  into  account  and  integra¬ 
tions  performed  for  more  appropriate  geometries,  detection  at  sensible  ranges 
might  then  be  predicted  by  this  theory.  No  effort  has  been  applied  to  this 
problem. 

3.2  Aerosol  effects  on  ambient  ionisation 

Rayleigh  scattering  of  radio  waves  from  the  volcanic  aerosol,  as  discussed  in 
Section  3.1,  is  the  only  available  mechanism  for  skywave  radar  detection  of 
dust  clouds  at  the  heights  where  they  could  pose  a  threat  to  aircraft.  The 
mechanisms  postulated  in  this  section  relate  to  phenomena  which  may  occur  at 
higher  altitudes  -  over  50  or  60  km,  say  -  and  in  general  it  would  not  be 
possible  to  infer  reliably  the  presence  of  low  altitude  dust  clouds  from  the 
observation  of  high  altitude  effects. 

3.2.1  Enhanced  recombination 

The  enduring  volcanic  aerosol  distributions  are  concentrated  in  the  lower 
stratosphere,  well  below  the  levels  of  strong  ionisation  (figure  12(ref.ll)). 
Accordingly,  the  most  likely  mechanism  for  influencing  radio  wave  propaga¬ 
tion  is  modification  of  D- layer  absorption  via  changes  to  the  electron 
density. 

At  D-region  heights  the  availability  of  three-body  channels 


X+  +  e  +  M  +X+M 


(15) 


and  the  presence  of  molecular  ions  with  the  degrees  of  freedom  to  undergo 
dissociative  recombination 


s 


.  »- 

SCU&S&i 
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XY  +  e  +  X  +'  Y 


mean  that  electron  loss  mechanisms  are  not  in  short  supply.  Below  80  km 
altitude  there  are  also  multiply-hydrated  proton  complexes  with  a  very 

high  recombination  coefficient  approaching  10-s  cm3  s'1 ,  compared  with 

...  M  ^ 

about  6  x  10  cm  s  for  NO  ions  which  dominate  above  80  km 
(figure  13 (ref. 12) ) . 

It  follows  that  for  aerosols  to  have  a  detectable  effect  some  particularly 
efficient  reaction  would  need  to  operate. 

The  reaction  postulated  here  is  catalytic  recombination  in  which  the 
surface  of  the  aerosol  particle  provides  a  100%  efficient  agency  for  the 
reactions 


X+  +  A  ♦  [  X+A] 


[  X+A]  +  e  *  [  XA]  *  X  +  A 


e  +  Z  +  Z" 


X  +  Z  ♦  X  +  Z 


In  order  to  estimate  the  rate  constant  for  this  process,  assume  the 
capture  cross-section  of  the  aerosol  is  equal  to  its  physical  cross- 
section.  Since  Tg  >  (figure  14)  the  inequality  vg  >?>  v^  is  ensured  and 

hence  the  rate  is  governed  by  v^.  From  simple  kinetic  theory  considera¬ 
tions,  for  diatomic  ions  which  dominate  below  100  km  (figure  15),  the 
recombination  coefficient  is  given  by 


a*  =  w 


5  kT. 

l 


N.  N  .  a 
i  a 


where  N.  and  N  are  the  ion  ar.d  aerosol  number  densities, 
x  a 

Now  T.  ~  T  .  ,  ~  3  x  102  °K  while  a  <  0.5  /am  above  40  km.  The  dominant 

x  neutral 

ion  species  is  NO+  so  nu  =  30  ra^;  N^  ~  10s  cm"3  at  90  km  (figure  15)  whence 
at  most. 


a*  ~  2.9  x  10"5  N  cm3  s"1 


V-f 
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Thus,  in  order  to  be  significant  in  competition  with  the  dominant  loss 


mechanisms  below  80  km  an  aerosol  particle  density  of  over  1  cm"3  is 
required. 

For  altitudes  in  the  range  20  to  40  km  where  aerosol  number  densities  near 

1  cm”3  could  apply  for  0.5  Mm  particles,  is  at  least  two  orders  of 

magnitude  less  than  it  is  at  90  km  (figure  16) .  On  the  other  hand,  for 

altitudes  near  the  D-E  transition  where  >  10s  cm"3 ,  aerosol  densities 

of  1  particle  (of  0.5  Mm  radius)  per  cubic  centimetre  equate  to  a  mass 

loading  of  10"6  g  cm"3  compared  with  an  air  density  of  about  10”7  g  cm"3 

(65  km)  or  10"8g  cm'3  (83  km)  (figure  11).  In  other  words,  the  air  is  not 
dense  enough  to  support  the  aerosol. 

Since  this  hypothetical  100%  efficient  recombination  mechanism  could  not 
produce  observable  effects,  it  is  reasonable  to  assert  that  more  realistic 
(and  more  complicated)  mechanisms  are  unlikely  to  do  so. 

3.2.2  Capture  of  free  electrons  without  recombination 

Another  hypothetical  process  which  could  be  advanced  as  a  possible  agent 
for  observable  effects  is  preferential  attachment  whereby  electrons  are 
readily  captured  by  aerosol  particles  according  to  equation  (19)  but,  for 
some  reason,  recombination  processes  such  as  equation  (20)  are  suppressed. 

In  the  simple  model  proposed  here,  the  aerosol  particle  accretes  charge 
until  electrostatic  repulsion  prevents  further  capture.  This  limit  corres¬ 
ponds  to  a  balance  between  the  thermal  energy  of  the  electrons  and  the 
Conlomb  potential  on  the  aerosol  surface: 


3  kT  Ne2 

2  ~  4xre  a 

o 


(23) 


At  altitudes  below  50  km  the  electron  density  is  too  low  to  allow 
appreciable  absorption.  For  the  upper  D- region,  Tg  ~  103  °K.  Then  the 

maximum  number  of  charges  which  could  be  accumulated  by  a  single  particle 
is  given  by 


N  ~  9  x  107  a 


(24) 


where  a  is  in  metres. 

Setting  a  =  0.2  Mm  which,  from  figure  8  seems  a  reasonable  upper  bound  for 
aerosols  above  50  km,  leads  to  N  =  18  electrons  per  particle.  Now,  for 


mid-D  region  and  night-time  E,  electron  densities  of  ~  103  cm'3  are 

typical.  It  follows  that  if  the  aerosol  density  were  to  approach  ~  102  cm 
a  large  fraction  of  the  free  electrons  could  be  captured,  thereby  removing 
the  main  absorption  mechanism  without  affecting  the  ionisation  density. 
Such  a  density  is  not  possible  at  a  =  0.2  Mm  because  the  corresponding 
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mass  loading  is  6  x  10'8  g  cm“3  which  is  ~  0.1  p(air)  at  50  km  but 
exceeds  p(air)  at  70  km. 

If  a  were  smaller,  say  0.02  pm,  the  charge  captured  would  decrease  by  one 
order  of  magnitude  but  there  would  be  103  more  particles  for  the  same 
density.  Using  a  mass  loading  of  0.1%  of  the  air  density  (equal  to  the 
mean  value  in  the  model  of  Section  3.1  based  on  Mount  St  Helens  data)  still 
leaves  the  depletion  a  factor  of  10  too  small  to  be  significant  at  70  km. 
Thus  only  a  minor  effect  on  absorption  appears  possible  even  when 
recombination  processes  are  prohibited. 

To  place  these  absorption  calculations  further  in  perspective,  pathloss 
experiments  during  the  JINDALEE  program  indicate  typical  two-way  absorption 
losses  of  about  6  dB(ref,13).  Thus  the  most  severe  effect  of  a  major 
reduction  in  non-deviative  absorption  is  still  essentially  a  threshold 
effect  given  the  other  variabilities  which  characterise  skywave  propagation. 

3.2.3  Depletion  of  sporadic-E 

There  is  no  available  data  on  aerosol  densities  at  E-region  heights  but 
any  reasonable  extrapolation  to  90  km  and  above  yields  exceedingly  low 
concentrations  of  only  the  smallest  particles.  Thus  ionisation  modifica¬ 
tion  of  normal  E-region  plasma  properties  can  safely  be  regarded  as 
negligible.  The  only  avenue  which  could  conceivably  lead  to  aerosol 
influence  is  interaction  with  sporadic-E  (Es) . 

Sporadic-E  is  believed  to  be  due  to  the  presence  of  long-lived  metallic 
ions  of  mainly  meteoric  origin.  The  principle  loss  processes  for  these 

ions  -  mainly  Mg  ,  Na  ,  Fe  and  Ca+  -  is  dissociative  recombination  via 
ion-exchange  reactions,  eg 


Mg+  +  02  *  Mg0+  +  0 


(25) 


Mg0+ 


+  e  * 


Mg  +  0 


(26) 


The  rate  coefficients  for  reactions  of  this  type  are  very  low, 

~  10  10  cm3  s  1 (ref . 14) .  It  could  be  the  case  that  active  sulphated 
radicals  produced  from  the  volcanic  aerosol  provide  a  depletion  mechanism 
for  these  metallic  ions  which  is  otherwise  lacking,  leading  to  the 
destruction  of  E^  layers.  From  the  point  of  view  of  detection,  however 

this  is  probably  irrelevant  since  the  stratospheric  aerosols  of  concern  to 
aircraft  would  not  be  correlated  with  the  high  altitude  phenomena  (eg 
figure  5).  In  any  event,  horizontal  structure  in  E  is  usually  so 

pronounced  that  observability  and  recognition  of  perturbed  E  would  not  be 
possible.  s 
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4.  CONCLUSIONS 

By  means  of  the  simple  physical  models  developed  in  this  paper,  rough  estimates 
have  been  derived  for  the  magnitudes  of  several  mechanisms  whereby  skywave  radar 
might  detect  volcanic  aerosols.  The  conclusion  in  each  case  is  that  detection 
is  not  viable. 

In  the  case  of  direct  radar  echoes  from  the  dust  clouds,  the  calculated  effects 
are  many  orders  of  magnitude  too  weak  to  be  detectable.  In  contrast,  some  of 
the  high  altitude  effects  appear  to  be  within  one  or  two  orders  of  magnitude  of 
being  observable,  but  these  effects  have  been  calculated  for  idealised  hypo¬ 
thetical  processes  which  may  not  occur  in  reality. 

Microwave  radar  detectability  of  dust  clouds  could  be  computed  by  a  simple 
extension  of  the  model  developed  here  for  HF  but  this  has  not  been  attempted. 


Postscript 

On  3  February  1984  it  was  learned  that  Mt  Galunggung  was  erupting.  The  JINDALEE 
radar  was  used  in  a  variety  of  modes  to  observe  the  region  around  the  volcano 
but  no  unusual  echoes  or  fading  could  be  detected. 
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TABLE  1.  COMPARISON  OF  18  MAY  ERUPTION  OF  MOUNT  ST  HELENS 
WITH  OTHER  ERUPTIONS (REF. 1) 


Mount 
St  Helens 


Worldwide 
maximum  past 
250  years 


Worldwide 
maximum 
all  time 


Total  magma  volume  erupted, 
(expressed  as  cubic  kilometres 
of  dense  rock) 

0.5 

>  1000 

,  '  .  '  . 

Eruption  magnitudea 

6 

10 

•  V:V->S' 

Eruption  column  height, 

23  to  27 

d 

a 

(.'iv:/::- 

(kilometres  above  volcano) 

maximum  18 

28  to  35° 

40  +  d 

■gwastfcafacEca 

Duration  of  climax  of  eruption 
(hours) 

average 

9 

38 

48 

(id 
*  .*  /, 

Volcanic  explosivity  index 
(VEI)b 

5 

7 

8 

. 

Area  covered  by  isopach 
enclosing  half  volume  of 
deposit  (square  kilometres) 

1  x  10s 

6  x  104 

106 

Area  of  complete  devastation 
(square  kilometres) 

500 

1000 

5  x  104  C 

’•,***  “•  * » 

a  Tsuya  (1955). 
b  Newhall  and  Self  (1982). 

C  Walker  (1980) . 

^  >  40  km  supposedly  reached  by  Krakatoa,  1883  and  Bezyraianny,  1956:  bot 
inadequately  substantiated. 
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TABLE  2.  ESTIMATES  OF  VOLATILES  RELEASED  BY  THE  18  MAY  ERUPTION  (REF.  1) 


Mass  of  magma  erupted  at  Mount  St  Helens 

18  May  volume3  0.3  km3 

post-18  May  volume  0.03  km3 


Total 


0.33  km3 


Initial  content 


<  500  ppm 
1000  ppm 
5  wt.  % 


Minimum  volatiles  erupted 

(0.33)  km3)  (250  ppm)  (2.4  g/cc)  x  101S  =  2  x  1011  g 
(0.33)  km3)  (1000  ppm)  (2.4  g/cc)  x  101S  =  0.8  x  1012  g 


(0.33)  km3)  (5%)  (2.4  g/cc)  x  101 


=  4  x  1014  g 


These  values  assume: 

(a)  no  meteoric,  hydrothermal,  or  atmospheric  water; 

(b)  no  scavenging; 

(c)  complete  degassing;  and 

(d)  no  contribution  from  intrusive  magma. 

Estimates  of  atmospheric  impact  will  have  to  consider  and  evaluate  these  factors 


Value  is  uncertain,  between  0.2  and  0.5  km3. 

Values  based  on  data  of  Melson  et  al  (1980). 
Note  particularly  the  uncertainty  for  S. 
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TABLE  3.  ERUPTION  RATES  AND  COLUMN  HEIGHTS  FOR  SOME  RECENT 
VOLCANIC  ERUPTIONS (REF. 1) 


Mount  St  Helens,  1980 
Hekla,  1947 
Agung,  1963 
Bezymianny,  19S6 
Santa  Maria,  1902 
Hekla,  1970 
Ngauruhoe,  1975 
Taupo ,  160  AD 


Maximum  volume 
flux  rate 
(m3/s  dense  rock) 


2  x  10’ 

2  x  104 

3  x  104 
2  x  10s 

4  x  104 
6  x  103 
2  x  103 
1  x  106 


Maximum  column  height 
(kilometres) 

above  volcanoa 


22(o) 
24(o), 21(c) 
23(c) 

45(o) ,42 (c) 
34(c), 29(o) 
16(o), 15(c) 
10(o), 10(c) 
>  50(c) 


o  =  observed,  c  =  calculated. 
b  A  vulcanian  (less  intense  than  plinian)  eruption. 

TABLE  4.  MAXIMUM  HEIGHTS  OF  ERUPTION  COLUMNS  AT  MOUNT  ST  HELENS (REF. 1) 

_  Maximum  column  heights  Duration  of  ash  emissic 

Date  of  eruption  (kilometres)  at  ahtwe  12 


Duration  of  ash  emission 
at  elevations  above  12  km 


18  May 

>  24.4,  17.3,  14.6,  19.2 

about  9  hours 

25  May 

12.2 

less  than  30  min 

12  June 

15.2,  10.7,  9.8,  10.7 

about  30  min 

22  July 

13.7,  14.5,  13.7 

25  min 

7  August 

13.4,  10.0,  6.1,  7.6, 

5.2 

2  to  5  rain 

16  October 

12.8 

less  than  5  min 

17  October 

14.3,  13.7 

less  than  5  min 

18  October 

5.2,  7.9 

less  than  5  rain 

«>£*.»  **7.1 


s-.--v-v.vi 

V; 

V-V-V-V- 
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TABLE  5.  VOLUMES  OF  SOME  RECENT  AIRFALL  DEPOSITS (REF. 1) 

Volume 

Vo  came  erup  ion  (cubic  kilometres  of  dense  rock) 


Mount  St  Helens,  1980 
Fuego,  1970 
Agung ,  1963 
Bezyroianny ,  1956 
Hekla,  1947 
Santa  Maria,  1902 
Krakatoa,  1883 
Tambora,  1815 
Taupo,  160  AD 


0.2  to  0.48 

0.1 

0.6  to  1.2  (March  and  May  events) 

0.5  to  2.0 

0.3 

4  to  9 

5  to  10 

30  (very  approximate) 

6  to  10 


I 
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Figure  2.  Stratospheric  aerosol  extinction  rates  measured  by  SAGE (ref. 4) 


CONCENTRATION  <cm'J> 


Figure  4.  Profiles  of  aerosol  number  concentration (ref. 3) 
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Figure  7 


I  I  18  TO  '9  km 

19  TO  20  km 

20  TO  21  km 


(a)  Map  showing  areas  and  altitudes  of  enhanced  extinction 
(50  percent  or  more  above  normal)  as  observed  by  SAGE 
between  April  21  and  28,  1979. 


-  FLIGHT  PATH 

CZZ)  AEROSOL  PLUME  BELOW  18  km 
Cm  AEROSOL  PLUME  AT  18  TO  19  km 
AEROSOL  PLUM  E  AT  1 9  TO  20  k  m 
I  I  NO  AEROSOL  PLUME 


<b)  Airborne  lidar  observations  in  the  vicinity  of  St.  Vincent, 
April  18  and  19,  1979.  Also  shown  are  the  calculated  air 
movements  at  three  pressure  levels  during  the  interval  between 
the  eruption  at  2057  GMT  on  April  17,  1979,  and  the  time  of  the 
flight. 


Figure  7.  SAGE  and  Lidar  measurements  of  Soufriere  Plumes (ref . 4) 


RELATIVE  MASS  CONCENTRATION  PER 
UNIT  LOGARITHMIC  DIAMETER  INTERVAL 
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Figure  11 


Figure  11.  Temperature  -  height  profiles  of  the  US  Standard  Atmosphere  1962 (ref. 10) 
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Figure  13 
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Figure  13,  Variation  of  diatomic  ion  concentrations  relative  to  0+(ref.l2) 
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Figure  14 


Figure  14.  Daytime  electron  and  ion  temperature  profiles (ref. 10) 
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